Equivalence of glass transition and colloidal glass transition in the hard-sphere limit 
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We show that the slowing of the dynamics in simulations of several model glass-forming liquids 
is equivalent to the hard-sphere glass transition in the low-pressure limit. In this limit, we find 
universal behavior of the relaxation time by collapsing molecular-dynamics data for all systems 
studied onto a single curve as a function of T/p, the ratio of the temperature to the pressure. 
At higher pressures, there are deviations from this universal behavior that depend on the inter- 
particle potential, implying that additional physical processes must enter into the dynamics of 
glass- formation. 
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In many amorphous systems, particle dynamics can 
slow down so dramatically as a control parameter is var- 
ied that fluid behavior is suppressed and the system forms 
an apparently rigid solid. One common example is a 
molecular liquid which, upon sufficient supercooling, be- 
comes a glass with no ability to flow Here, the tem- 
perature T governs, among other things, what energy 
barriers may be crossed as the system attempts to equi- 
librate @. Another example is a colloidal fluid, in which 
micron-sized particles driven by Brownian motion form 
a colloidalglass when the density </>, or pressure p, is 
raised 0, Hf. In this case, density controls the parti- 
cle dynamics by governing the amount of free volume 
available for particles to rearrange While the con- 
ceptual frameworks for understanding the slow dynam- 
ics in these two cases appear to focus on two very dif- 
ferent physical effects, many aspects of these structural 
glass transitions appear to be similar Q . In neither case 
is it known whether an underlying phase transition ac- 
companies glass formation; as the dynamics slow down, 
no obvious structural changes or unambiguous static or- 
der parameters emerge [H, H|. Without such a theoret- 
ical underpinning, it has been impossible to determine 
whether both phenomena are driven by the same physics 
or whether varying T and p slows the dynamics in intrin- 
sically different ways. 

By analyzing relaxation in several models that have 
frequently been used to capture different features of dy- 
namical slowing, we are able to identify a universal aspect 
of glass formation that emerges at low pressures and tem- 
peratures. In this limit, we can scale the relaxation times 
r for all systems studied onto a single master plot as a 
function of the single variable T/p. This has two impor- 
tant consequences. On the one hand, the scaling collapse 
allows one to see unambiguously that at sufficiently low 
pressures and temperatures, the colloidal glass transition 
and molecular glass transition are manifestations of the 
same phenomenon. In this limit the temperature-driven 
molecular glass transition is equivalent to the pressure- 
driven hard-sphere colloidal glass transition; temperature 



and pressure are equally important for slowing relaxation 
processes. On the other hand, deviations from this uni- 
versal behavior away from this limit demonstrate that at 
least two distinct physical processes can enter to produce 
the slowdown of particle motion. 

Our models consist of spheres with interparticle inter- 
actions V(rij), where is the distance between particles 
i and j with diameters Oi and cr, , respectively. The par- 
ticles do not interact at large separations: V(ry) = 
for Tij > (Tij, where = (erj + <Jj)/2. If r ?J < cry, the 
particles repel each other according to either the Weeks- 
Chandler- Andersen potential @ 



v( rij) = - 

or a simple power law 



-(1-r^/a^r, 



(1) 



(2) 



where we use a = 2 (harmonic repulsion), a = 5/2 
(Hertzian repulsion) or a = (hard-sphere repulsion). 
All of these systems have a jamming transition in the 
limit of zero temperature, T — ► 0, and zero pressure, 

We use molecular dynamics at fixed temperature and 
pressure to simulate three-dimensional systems composed 
of 50:50 mixtures of N = 1000 particles with diameters 
a and ar, = lAa and the same mass, m. Temperature 
has units e, pressure has units e/cr 3 , and time has units 
\Jm<j 2 / e. The Boltzmann constant, fcs, is set to unity. 
We use an event-driven code for the hard-sphere simu- 
lations. Periodic boundary conditions arc applied in all 
directions. We calculate the self-part of the intermediate 
structure factor: S{k,t) = -j| J2i exp(zfc-[rj(£)— ^(0)]) [1], 
where the sum is over all large particles, ri(t) is the lo- 
cation of particle i at time t, and k is chosen in the 
x— direction. The amplitude of k satisfies the periodic 
boundary conditions and is approximately the value at 
the first peak of the static structure factor. We define 
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FIG. 1: Relaxation time versus different control parameters 
for a system of 1000 particles with harmonic repulsions, (a) 
Relaxation time r versus temperature T at different fixed 
pressures, p : p = 2 x 10~ 7 (circles), p — 2 x 10 -6 (squares), 
p = 2 x 10~ 5 (diamonds), and p = 2 x 10~ 4 (upward tri- 
angles), (b) Relaxation time r versus inverse pressure \jp at 
different fixed temperatures: T = 10 -8 (downward triangles), 
T = 10~ 7 (pluses), T = 10~ 6 (crosses), and T = 10~ 5 (stars). 



the relaxation time r to be the time at which t) = 
e^ 1 S(k, 0). We take data after a system has been equili- 
brated for several r. 

Fig. [TJa) shows the relaxation time r versus temper- 
ature T for a system with harmonic repulsions (a = 2 
in Eq. [2]) that is cooled at different fixed low pressures, 
p. This corresponds to the standard trajectory for ex- 
periments on supercooled liquids undergoing the glass 
transition, a trajectory in which p is typically fixed at 
atmospheric pressure. Fig. QJb) shows t versus 1/p at 
different fixed T for the same system; in these trajec- 
tories, wc raise p at fixed T as is typically done in ex- 
periments on colloidal systems. As expected, r increases 
with decreasing T and increasing p. 

Fig. HJa) shows that we can collapse the data in Fig. [1] 
for all the trajectories, both at fixed pressure and fixed 
temperatu re, ont o a single curve by scaling the relaxation 
time by ^Jm/pa and temperature by pa 3 , so that 



-/sjm/pa = F(T/pa 3 



(3) 



This collapses data ranging over 4 decades of temperature 
and pressure. This remarkable collapse is achieved by 
using the time scale, ^Jm/pa, and energy scale, pa 3 , to 
make relaxation time and temperature dimensionless Q . 
The characteristic time \J m /pa is proportional to the 
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FIG. 2: Collapse of all the relaxation time data shown in 
Fig- H ( a ) Scaled relaxation time, r/y^m/pa, versus scaled 
temperature T/pa 3 . (b) Scaled relaxation time r/ ym /pa 
versus scaled pressure pa 3 /e, for T/pa 3 = 0.08 (circles), 
T/pa 3 = 0.1 (squares) and T/pa 3 = 0.2 (diamonds). Hor- 
izontal lines show the limiting values of r^/pa/m as pa 3 /e — > 
0. 



time for a particle starting at rest to move its diam- 
eter a due to a pressure p and is the duration of a 
pressure-driven particle rearrangement. The dependence 
on T/pa 3 implies that both T and p are equally impor- 
tant in controlling the dimensionless relaxation time. 

Dimensional analysis provides a starting point for un- 
derstanding the implications of the d ata collapse of Eq. [3] 
The dimensionless time r / ^Jm/pa can be written as a 
function of the dimensionless variables of the system. In 
addition to T/pa 3 , another dimensionless ratio is pa 3 /e, 
where e sets the scale of the interaction energy in Eq. [2] 
There are no other independent dimensionless variables 
for the system with harmonic repulsions. The dimension- 
loss relaxation time must therefore satisfy 



-^pa/m = f(T/pa 3 ,pa 3 /e 



(4) 



The data shown in Fig. [UJa) all lie at low pressures, 
where the second argument in Eq. [51 pa 3 /e, is small. 
In Fig. [2£b), we show the scaled relaxation time versus 
pa 3 /e, at three values of T/pa 3 . In all cases, the data 
approach an asymptotic value at low pa 3 /e. Thus, in the 
low-pa 3 / e limit, r^Jpa /m is a function of T/pa 3 only, 
consistent with the collapse of Fig. [2ja). 

The limit, pa 3 /e — > 0, is always satisfied in the hard- 
sphere limit, e — > oo. Thus, the relaxation time for 
hard spheres should collapse onto the same scaling form 
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FIG. 3: (Color online) Scaled relaxation time, t/ \Jm/p<r, 
versus scaled temperature T/pa 3 for all the data for the har- 
monic potential (black) in Fig.[T]as well as for the Hertzian po- 
tential, a = 5/2 (red), the hard-sphere potential (magenta), 
and the Weeks-Chandler- Andersen potential (blue). Black 
solid curve is the Vogel-Fulcher fit: y = 0.59 exp(0.I8/(a; — 
0.045)), where y — r/ ^Jm/pa and x = T/pa 3 . Blue- 
dashed curve is a fit to the Elmatad-Chandler-Garrahan form: 
y = 3.1 exp[0.064(a:- 1 - 4.72) 2 ]. 



as in Fig. [5] Moreover, all potentials that behave as 
hard spheres in the low-pressure limit by preventing over- 
lap between particles — namely, all potentials with finite- 
ranged repulsions-should also collapse onto the same 
form at low pa 3 /e. This is not what we would expect 
from dimensional analysis alone, since different poten- 
tials can contain additional dimensionless parameters, 
such as the exponent a in Eq.[^l However, the physics of 
the pa 3 /e — >• 0, hard-sphere limit suggests that collapse 
should occur, irrespective of the value of a, as long as 
a > so that the potential is repulsive. 

This is corroborated in Fig. [3J where we show the data 
of Figs. [T] and [2] for harmonic repulsions (a — 2 in Eq. [2]) 
together with data for three other potentials: Hertzian 
repulsions (a = 5/2 in Eq. [5]), the hard-sphere poten- 
tial (a = in Eq. [5]), and the Weeks-Chandler- Andersen 
potential (Eq. [T]) . Indeed, all of these systems collapse 
onto the same scaling form. Thus, for the low-pressure 
(low pa 3 /e) data considered here, r^/pa /m is a function 
of T/pa 3 only and does not depend separately on the 
interaction potential - cither on its form (i.e., whether 
it is given by Eq. [T] or [2] with different exponents a) 
or on its overall magnitude, e. This represents a major 
simplification of the relaxation-time data. 

While the scaled dynamics are independent of dimen- 
sionless parameters characterizing the inter-particle po- 
tential, they may be affected by other dimensionless num- 
bers characterizing the system. For example, relaxation 
time data for hard-sphere systems with different polydis- 
persities or diameter ratios will not necessarily collapse. 

The data collapse suggests a way of looking at the glass 
transition in this limit where the relaxation time is only 
a function of T/pa 3 . The ratio of T/p corresponds to an 
effective volume created by using thermal energy to do 
work against the pressure: pAV ~ T. In the hard-sphere 




FIG. 4: (Color Online) Equation of state for three- 
dimensional systems at low pa 3 /e for hard sphere systems 
(solid magenta circles) and systems with harmonic repulsions 
measured along both constant temperature and constant pres- 
sure trajectories. The dashed line is a fit to free- volume the- 
ory, pa 3 /T = 0.980/(1 — (0/</> c ) 1/3 ), with <j> c = 0.66. The solid 
line is an empirical fit given by pa 3 /T = 0/(1 - (0/1. 13)) 4 35 . 



glass transition, the system relaxes via free volume. Like- 
wise, soft-sphere liquids can create free volume by using 
thermal energy. At high temperatures, where there is 
plenty of this "thermal" free volume, the system relaxes 
rapidly; at low temperatures, where there is less thermal 
free volume, the system relaxes more slowly. 

The form of the scaling function, F(x), shown in Fig.[3J 
should tell us whether the system has a thermodynamic 
glass transition. If such a transition exists, as has been 
suggested theoretically for hard spheres [l(|, then F(x) 
should diverge at a nonzero value of T/pa 3 , i.e., at 
nonzero T for soft-sphere liquids or at finite p for hard 
spheres. Fig. [3] shows that over the dynamic range of our 
simulations, the Vogel-Fulcher form, 



r / y m/pa = C exp(A/ (x — xq)) 



(5) 



with x = T/pa 3 , provides a reasonable fit (solid line) 
with x = 0.045, A = 0.18 and C = 0.59. The scaling 
collapse and Eq. [5] imply that To /pc 3 = xq so that To = 
pa 3 xo- In other words, To increases with pressure. These 
results are consistent with recent numerical studies on 
spheres with harmonic repulsions [Til ] . which show that 
To increases with packing fraction above some critical 
value. The parameter Tq is often used to parameterize 
the fragility of the system [l^, HH; thus, the fit to Eq. [5] 
implies that fragility increases with pressure. 

We also find that a fitting form proposed by Elmatad, 
et al. [3], for which r/ y/m/pa diverges only at T/pa 3 = 
0, provides an equally good fit in the regime of interest 
at small T/pa 3 . For this fitting form, 



<- / \J m/pa = Ci exp 



1 



X\' 



(G) 



where x = T/pa 3 . We find Xi = 0.21, A 1 = 0.25 and 
C\ =3.1. Thus, we cannot distinguish whether the scal- 
ing function diverges at nonzero or zero T/pa 3 . This is 
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not surprising, since even experiments with 17 decades of 
dynamic range cannot tell whether there is a thermody- 
namic glass transition. However, we note that if there is 
no thermodynamic glass transition, the scaling function 
diverges at T/pa 3 = 0, so that the relaxation time di- 
verges in the double limit T/pa 3 — ► 0, pa 3 /e — ► 0. This 
double limit corresponds to the zero-tempcrature jam- 
ming transition of frictionlcss spheres with finite-ranged 
repulsions, also known as Point J 0. This implies that 
the glass transition is controlled by Point J if there is no 
intervening thermodynamic glass transition. 

The relaxation time is not the only quantity to exhibit 
data collapse in the low pa 3 /e limit. Other unrelated 
quantities, including those that are independent of the 
dynamics, such as the packing fraction, <j>, should also 
exhibit collapse in that limit. The packing fraction is the 
number density made dimensionless by the average par- 
ticle volume. For a fixed potential, <f) must be expressible 
as a function of T/pa 3 and pa 3 /e. In the low-pressure, 
hard-sphere limit where pa 3 /e < 1, it should satisfy 
qt> = H(T/pa 3 ) with the same scaling function H(x) for 
all finite-ranged repulsive potentials. This function can 
be inverted to yield the equation of state 

pa 3 /T = H (</>). (7) 

Fig. 2] shows the data collapse for the equation of state 
along different trajectories and for different potentials for 
our bi-disperse systems in three dimensions. We have 
shown two fits to the data. The dashed line is a fit to 
free volume theory for (j) near the fitting parameter <j) c — 
0.66. This form fits experimental data for colloidal hard 
spheres and numerical data for hard spheres reasonably 
well but is clearly unsatisfactory here. The solid line 
is an empirical fit to the data (see caption). 



Our results show that the data for the relaxation time 
in many different systems collapse onto a single curve of 
Tyjpa jra versus T/pa 3 . This collapse, however, is con- 
fined to p <C e/a 3 . We find that upon increasing the pres- 
sure, there are deviations from the scaling collapse (see 
Fig. EJb)). This implies that there are additional contri- 
butions to relaxation near the glass transition, beyond 
the particle rearrangements facilitated by the thermal 
free volume. At large pa 3 /e, relaxation also occurs be- 
cause thermal fluctuations drive the system across energy 
barriers. This contribution could perhaps be taken into 
account by introducing an effective pressure-dependent 
radius for the particles (such as the Barker-Henderson 
radius [IB])- As pa 3 /e increases at fixed T/pa 3 and the 
density of the system increases above the jamming point, 
the effective radius shrinks, increasing the effective free 
volume beyond the thermal free volume created by the 
thermal energy working against the pressure. 

The relative effects of T and p on r in glass-forming 
liquids have been studied experimentally, with differing 
conclusions EH ■ I n order to understand the conse- 
quences for molecular liquids, we must first understand 
the corrections to the leading hard-sphere behavior when 
pa 3 / e is no longer small. Molecular liquids typically have 
densities high compared to that of a system near the jam- 
ming threshold. They also have long-ranged attractions 
as well as the short-ranged repulsions considered here. 
For these glass-forming liquids, hard spheres may still be 
a useful starting point, but at least one other distinct 
contribution to the relaxation must also be considered. 
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